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Abstract 

We measure hole transport in poly(3-hexylthiophene) field effect transistors with channel lengths 
from 3 ^m down to 200 nm, from room temperature down to 10 K. Near room temperature effective 
mobilities inferred from linear regime transconductance are strongly dependent on temperature, 
gate voltage, and source-drain voltage. As T is reduced below 200 K and at high source-drain 
bias, we find transport becomes highly nonlinear and is very strongly modulated by the gate. We 
consider whether this nonlinear transport is contact limited or a bulk process by examining the 
length dependence of linear conduction to extract contact and channel contributions to the source- 
drain resistance. The results indicate that these devices are bulk-limited at room temperature, and 
remain so as the temperature is lowered. The nonlinear conduction is consistent with a model of 
Poole-Frenkel-like hopping mechanism in the space-charge limited current regime. Further analysis 
within this model reveals consistency with a strongly energy dependent density of (localized) valence 
band states, and a crossover from thermally activated to nonthermal hopping below 30 K. 

PACS numbers: 
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I. INTRODUCTION 



Electronic devices based on organic semiconductors, such as polymeric field effect tran- 
sistors (FETs) and light emitting diodes (LEDs), have attracted much interest as possible 

nn 

inexpensive and flexible alternatives to inorganic devicespj, While there has been con- 
siderable improvement in device properties, the detailed mechanism of electronic transport 
in organic thin film devices remains a subject of active investigation. 

Charge motion in undoped organic field effect devices is often characterized in the linear 
(small source-drain bias) regime by a field-effect mobility. In disordered polymer semiconduc- 
tors this mobility, fi, is strongly temperature dependent near room temperature, consistent 
with thermally assisted hopping between localized states dispersed throughout the polymer 
filmP]. These localized states are likely to be polaronici4,|- Several different models for 
hopping transport in these materials have been used to interpret experimental data from 
~200 K to room temperature, including simple thermal activation[^, 2d variable range hop- 
ping (VRH)j^, and percolative VRH with an . While some 
models work better with some specific samples, generally distinguishing between them is 
difficult. 

The situation is complicated by the fact that the effective mobilities inferred in this man- 

n 

ner depend on T, gated charge density |8|, source-drain bias, and contact effects. Parasitic 
contact resistances in particular can be important. Experiments in several organic semi- 
conductor/electrode combinations 0, [l^ 111, 12, [3, JJ, 13, ll^ have shown that contact 
resistances can be a significant fraction of the total source-drain resistance in the linear 
regime in short devices. Data in this shallow channel limit must be corrected accordingly 
for this parasitic series resistance, Rg, to find the true mobility within the semiconductor. As 
with /i, one must bear in mind that Rg typically depends on temperature, gate voltage, and 



local electric field 



12j |. When examining transport properties of organic FETs, it is therefore 



important to determine whether the devices are dominated by the bulk (channel) or the 
contacts. 

We report nonlinear transport measurements in field-effect devices made from high qual- 
ity, solution cast, regio-regular poly(3-hexylthiophene) (P3HT), with channel lengths from 
3 fim to 200 nm and aspect ratios of 10. Higher temperature properties at low source-drain 
fields are consistent with those observed by other investigators. From 200 K to 10 K, we ob- 
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serve gate-modulated nonlinear IV characteristics. To understand the effects of contacts in 
our series of devices, we examined both the devices described above, and an additional series 
of fixed-width FETs to obtain the channel and contact resistances as a function of temper- 
ature. These data demonstrate that (a) the fixed-apect-ratio devices are bulk (not contact) 
limited at high temperatures; and (b) the contribution of contacts relative to the channel 
actually decreases as the temperature is lowered, so that bulk-limited devices tend to remain 
so as T is decreased. Since the nonlinearities in the IV curves become more pronounced 
at low temperatures and in wider devices (for which contact resistance are proportionately 
less important), it is unlikely that these nonlinearities are due to contact effects in these 
geometries. We find that the nonlinear data are consistent with a model of Poole-Frenkel 
(PF) type conduction in the space-charge limited (SCL) regime. Within this framework, the 
strong gate and temperature dependence of this conduction are consistent with a density of 
localized valence states that varies approximately exponentially in energy. Still within this 
model, at temperatures below 30 K there appears to be a crossover from thermally assisted 
hopping to a nonthermal mechanism. These nonlinear data over a broad temperature range 
constrain any other models of transport in such devices. Studying FETs in this nonlinear 
regime allows comparisons between models not readily performed with linear transport data. 



II. EXPERIMENTAL DETAILS 

Devices are made in a bottom-contact configuration (see Fig. ^ inset) on a degenerately 
doped p+ silicon wafer to be used as a gate. The gate dielectric is 200 nm of thermal 
Si02. Source and drain electrodes are patterned using standard electron beam lithography. 
The electrodes are deposited by electron beam evaporation of 4 nm Ti and 25 nm of Au 
followed by liftoff. The fixed-aspect-ratio devices have channel lengths ranging from 3 fim 
down to 200 nm, with the channel width scaled to maintain w/L = 10. Larger FET devices 
{w = 1 mm, L = 50 fim) are also prepared for comparison. We also examine a second set of 
devices with fixed width w = 100 fim and channel lengths varying from 5 fim up to 40 /im, 
specifically for probing contact resistance issues. 

The organic semiconductor is 98% regio-regular P3HT0], a well studied material jisl liol 
I20I . P3HT is known from x-ray scattering to form nano crystalline domains with sizes on 
the order of 20 nmjisl. l2ol |. and more ordered films correlate with higher measured mobili- 
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ties. RR-P3HT is dissolved in chloroform at a 0.02% weierht concentration, and is solution 
oastQ o,.o cone-cleaned, eholo.ofo^.wabbed substrates^ The fixed-a.peet.atio se.e. 

of devices are from one casting, while the fixed-width series are from a second casting. The 
resulting film thicknesses over the channel region are tens of nm as determined by atomic 
force microscopy (AFM), though casting produces somewhat nonuniform films. All devices 
are stored in vacuum dessicators until use. The measurements are performed in vacuum 
(~ 10~® Torr) in a variable-temperature probe station using a semiconductor parameter 
analyzer (HP4145B). 



III. RESULTS AND DISCUSSION 



At room temperature, the devices operate as standard p-type FETs in accumulation 
mode liol 12^. 21 1. With the source electrode as ground, in the linear regime we 

extract p| an effective mobility from the transconductance. That is, from data of source- 
drain current. Id versus the gate voltage, Vg, at a fixed low drain voltage, Vd, we compute 
= (gmL) / {wCiVn) , where gm = dla/dVc is the transconductance, is capacitance per 
unit area of the gate insulator, and vol L is the aspect ratio. The relevance of parasitic 
contact resistances will be addressed below. As is reported elsewhere |y!], the mobility is 
gate voltage dependent, increasing with increasing Vq- It also increases with increasing 
source-drain voltage. Effective mobilities are typically between 10~^ and 10~^ cm^/Vs, and 
apparent threshold voltages (V^), though not necessarily meaningful j2l|. are low (< 2 V). 
For L = 50 /im FETs operated in the saturation regime, the on/off ratio is typically ~ 650, 
comparing between gate voltages of -95 V and V. As temperature is reduced below 300 K, 
the off-current drops to undetectable levels by 150 K, as the unintentional carriers (due to 
slight doping from air exposure) freeze out. 

Over the moderate temperature range of 300 K to ~ 200 K, the mobility as inferred 
above at fixed small Vd is found to depend steeply on temperature. A small representative 
set of this data is shown in Fig. ^ where /z is plotted vs. inverse temperature for Vq = 
—30, —40, and -50 V, for the L = 300 nm, w = 3 fim device at constant source-drain 
electric field of 1.3 x 10'' V/m. 

These data are approximately equally consistent with the three models mentioned above: 
simple Arrhenius behavior (the dashed line) with an activation energy ~ 100 meV; VRH 
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for a 2-D system (dotted line), of the form /i = /io exp(— (Tq/T)^/^), where /io and Tq 
are fit parameters; and finallvthe more sophisticated percolative VRH theory (sohd hnes) 
developed by Vissenberg et a/p,^. The Vissenberg model's underlying assumptions include 
an exponential density of (localized) states (DOLS), z/(e) ~ exp(e/A;BTo), with transport of 
carriers dictated by percolative hopping. Since the gate voltage controls the Fermi level in 

inds that transport in 



2l| are: Tq, describing 



the channel, and hence the occupation of the localized states, one 

u 

the channel is strongly affected by Vq- Relevant fit parameters j7|, 
the energy dependence of the DOLS; a prefactor ctq with units of conductivity; and a, an 
effective overlap parameter for hopping. Values used in the fits shown are Uo = 7 x 10^ S, 
To = 418 K, a = 4.35 x 10^ m~^, consistent with those seen by other investigators in 



P3HT 



As temperatures are lowered from 200 K down to 10 K, over a broad range of source- 
drain and gate voltages. Id evolves from approximately linear to a strongly nonlinear (su- 
perquadratic) dependence on Vd- An example of this evolution is shown in Fig. [21 for the 
L = 3 fim, w = 30 fim device, comparing data at 300 K and 70 K. We note that, at the low- 
est temperatures, smaller devices transport current more easily (larger currents at smaller 
gate voltages for a fixed Vd/ L) than larger devices, as we will discuss later. 

Analysis below shows that the nonlinear IV characteristics are decribed well by a model 
incorporating space-charge limited currents, modified by a Poole-Frenkel-like exponential 
dependence of effective mobility on square root of the local electric field (SCLPF). This 



conduction mechanism has been seen repeatedly in two-terminal devices 



perature experiments 
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23l | . Room tem- 



251 ] on P3HT FETs with 70 nm channel lengths also show indications of 



SCL currents. We find that within this model, the temperature and gate voltage dependence 
of the data support a strongly energy dependent DOLS such as that in the Vissenberg pic- 
ture. Other models may be possible, but they are constrained by the dependences presented 
below. 

Charge transport in a device is space charge limited if the injected carriers significantly 
alter the local electric field from the average field imposed by the electrode potentials, and 
correspondingly limit the current. If, instead, the bottleneck in charge transport is injection 
at the contacts, a device is said to be contact limited, and is expected to exhibit Ohmic 
behavior at low source-drain fields. In principle, modeling our devices requires the full 
solution of the steady state charge and electric field profile in a three-terminal accumulation 
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mode transistor, including field- and temperature dependent effective mobility, and field- 
and temperature dependent contact properties. This general problem is very complex j^]; 
here we consider a simpler model and compare with the transport data. 

In a system that is not contact limited, when the effective mobility varies as a function 
of electric field F, the space-charge limited current in a two-terminal device (a Id model) is 
determined by the numerical solution of 



J = p{x)efi{F{x))F{x) 



KeQ dF 

— = p[x) 

e ax 



(1) 



where p{x) is the local hole density, jj" F{x)dx = V, and the appropriate boundary condition 
on F{x = 0). For approximately Ohmic contacts, F{x = 0) ~ 0. 
We note that the case of an effective mobility that varies as 

/i(F) = /ioexp(7yF) (2) 

and Ohmic contacts has been solved j^l; and that the exact numerical solution is very well 
approximated by: 

Q /T/\ 2/1 r /T/\ 

(3) 



^exp 



V 

0.97 ( J 



Here A is the device cross-sectional area, and L is the interelectrode distance. In our geome- 
try L is the channel length, A is an effective cross-sectional area for the device (proportional 
to channel width), and k is the relative dielectic constant of the semiconductor (chosen to be 
3 in our analysis). The appearance of V/L in this equation does not imply that the electric 
field is constant over the device length. Rather, Eq. El suggests a means of plotting IV data 
to quickly ascertain consistency with the detailed numerical solution to the two-terminal 
SCLPF problem. 

Note that /xq can depend on temperature; in a picture of hopping it should be proportional 
to the effective DOLS at the injecting contact. The dependence of mobility on exp(\/F) 
has long been seen in semiconducting and conducting pol ymers |28l| . and is associated with 
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29, boll. The numerical solution 



t.ec.a.^c™.a.d diso.de. in these ,nate«sB 
of Ref. |27| (approximated by Eq. should be valid as long as the functional form of the 
field-dependent mobility remains exp(V^). 



This equation is derived |23l. I^^l assuming that the charge distribution is determined by 



source-drain electrostatics only. The electric field from the gate certainly plays a nontrivial 



role in our devices, clearly affecting charge injection, and allowing the formation of a channel 
at higher temperatures (since any charge present in the channel below dopant freeze-out has 
to be injected from the source and drain). For simplicity, however, in this model we will 
assume that the gate dependence will manifest itself through fiQ, and that for fixed gate 
voltage we may treat the source-drain conductance like a two-terminal device. 

FigureEt' shows a representative log-log plot of Id vs. Vd for one sample with L = 500 nm 
and Vg = —75 V, for several temperatures. The solid lines are the numerical solution to 
Eq. HJ with parameters A^q and 7 chosen at each temperature to give the best fit. The 
numerical solution is virtually indistinguishable from the analytic form of Eq. IHl Data on 
this and other samples for different gate voltages are qualitatively similar, with very good 
agreement between the numerical solution of the SCLPF model and the data. Fig.lSb shows 
the same data and fits replotted as suggested by Eq. Q. The quality of this agreement 
between the SCLPF model and the data in several devices over a broad range of T, Vq, and 
Vd/L is striking. Clearly for a given (large) Vg and (low) T, Id ~ exp(-\/V^). Alternate 
explanations of the nonlinear conduction are strongly constrained by this dependence. The 
temperature and gate voltage range over which this form of source-drain nonlinearity occurs 
varies systematically with sample size, as described below. 



A. Contact effects 



One must consider whether the nonlinear Id — Vd characteristics result from nonlinear 
contact resistances as the temperature is decreased. For several reasons, discussed below, 
we do not believe this to be the case. 

First we consider directly inferring the contact and channel resistances {Rs and -Rch, 
respectively) in the linear regime, and examining the temperature variation of their relative 
contributions. For a series of devices with fixed width, these resistances are calculated as 



follows 



|15|,Q: the total resistance, -Ron = ^Id/OVd, is calculated for 

each device at a small Vd and is plotted as a function of L for each gate voltage. The channel 
resistance per unit length, Rch/L, at a given T and Vg is the slope (dRon/dL) of such a 
graph, and the intercept {Ron extrapolated to L = 0) gives the parasitic series resistance, 
Rs. 

For a fixed- aspect-ratio series of devices, one may follow an analogous procedure. The 
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total source-drain resistance Ron = {L/w)Rci + Rs, where i?n is the resistance per square of 
the channel At a given T and Vq, Ron x w is plotted versus L for the series of devices. The 
slope of such a graph gives R\j, and the intercept gives the total parasitic contact resistivity, 
Rs X w. This analysis is shown in the fixed-aspect-ratio devices in Fig.EJat room temperature 
for several gate voltages. For our geometry of w/L = 10, the inset shows the inferred Rs/Rch 
as a function of Vq, for a device with L = 1 fim. Our series of fixed- aspect-ratio devices is 
clearly channel-limited at room temperature. 

This analysis may be repeated at different temperatures to examine the evolution of 
Rs and Rch- We find that the nonlinear conduction at large average source-drain fields 
shown in the previous section makes this difficult to measure over a broad temperature 
range in the fixed-aspect-ratio device series. However, the fixed-width devices with longer 
channel lengths are well-suited to this approach down to 100 K. Figure El shows Rs/Rch as 
a function of temperature for the L = 5 fim device from the w = 100 fim series. Near room 
temperature, -Rs > -Rch for this device. As T decreases, while both Rs and -Rch increase 
significantly, Rs falls below -Rch near 100 K. The results of this experiment and others to be 
published elsewhere demonstrate that the contact contributions become less important 
at low temperatures. This strongly suggests that the nonlinearities in the Id — Vd curves 
of short channel length devices observed at lower temperatures are unlikely to be contact 
effects. 

Furthermore, the trends in transport with sample dimensions also support this conclusion. 
We find at low temperatures that the smallest devices actually transport charge considerably 
better than larger devices. For example, at 50 K for fixed Vd/L, the fixed-aspect-ratio 
L = 300 nm device exhibits measurable conduction for gate voltages as small as -25 V, 
while the L = 3 fim device requires Vg = —45 V. This trend is the opposite of what one 
would expect for contact-limited conduction [ll|. Since w/L is held constant in this set 
of devices, shorter channel devices have significantly smaller contact well, further 

emphasizing this point. The data are, however, consistent with the suggestion of space charge 
effects seen in 70 nm channel length P3HT FETsQ- 

Finally, the voltage and temperature 
dependence of the data in Fig. El is not consistent with the forms for either classical Schottky 
contacts or Fowler- Nor dheim emission. Detailed theory ^ and experiments 32] on injection 



into disordered polymer semiconductors show that injection efficiency can actually improve 
as temperature is decreased, consistent with our contact resistance data described above. 
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Coupled with the size dependence, this supports the idea that low temperature transport in 
our devices is bulk limited rather than contact limited. 



B. Physical significance of fit parameters 

Continuing within the SCLPF model and our analysis of the fixed-aspect-ratio devices, we 
note that, for identical effective mobilities and 7 parameters, Eq. IHlimplies that two samples 
with the same aspect ratio, thickness, and average source-drain field should give the same 
currents, independent of channel length, even deep in the nonlinear regime. Device-to-device 
variability in the effective mobility and 7, presumably due to differences in P3HT thickness 
and microstructure, make this challenging to check directly in our devices. For reasonable 
values of 7 and Vd/L, a 10% variation in 7 would lead to more than a factor of 2 variation 
in predicted current at low temperatures because of the exponential dependence in Eq. 01 
However, if one fixes Vd/L and Vq-, and normalizes measured currents by room temperature 
mobilities, one does indeed find scaling. For example, the L = 500 nm and L = 1 /xm 
currents in the nonlinear regime normalized this way agree well all the way down to 10 K. 

Within this model, the parameter 7 should depend only on the hopping mechanism {e.g. 
thermal activation) and the nature of the localized states. We therefore expect 7 to be 
independent of gate voltage for a given sample, and this is indeed seen in the inset to 
Fig.ini At high temperatures (T > 50 K), the data for all gate voltages and all samples look 
roughly linear in 1/T, consistent with thermally activated hopping. The magnitude of the 
slope of 7 vs. 1/T is approximately 0.12 K(m/V)^/^. This is consistent in magnitude with 
coefficients found in other semiconducting polymers such as poly(phenylene vinylene) 22 1. 
However, 7 vs. 1/T deviates significantly from a straight line at lower temperatures for all 
samples. This is consistent with a crossover from thermally activated hopping to a much 
less steep temperature dependence. A natural candidate is field enhanced tunneling between 
the localized states. 

We now consider the gate and temperature dependence of the parameter AfiQ found by 
the numerical analysis above. The effective cross-section for current fiow. A, is assumed 
to be temperature and gate voltage independent for each sample. In the linear regime at 
moderate temperatures it is known that the mobility inferred from the transconductance is 
gate voltage dependent, as seen in Fig. [H This dependence on Vq continues in the apparent 
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SCLPF regime, as shown in Fig. [7^ on an L = 300 nm device. The amount of variation of 
fiQ with Vq, roughly a 5% increase of fio per voh of Vg for this sample, shows no strong trend 
with temperature. The magnitude of this variation of mobility with Vq is consistent with 
that seen at higher temperatures in Fig. Q for this sample. This exponential dependence 
of fiQ on Vg is seen throughout the apparent SCLPF regime. As in the linear regime of an 
accumulation FET, gate voltage modulation of the conduction is unipolar, higher currents 
result only when Vq is made more negative. 

The temperature dependence of /io is also very strong, as shown in Fig. [7|d for the same 
300 nm channel sample at five gate voltages. The variation of fJ^oiT) shown is much closer 
to an exponential in T to some slightly sublinear power rather than an Arrhenius or VRH 
form. This strong temperature dependence of fiQ holds for all samples in the apparent 
SCLPF regime. 

The strongly energy dependent DOLS employed in, for example, the Vissenberg model 
offers a natural explanation for these steep dependences of /iq on Vg and T. In the absence 
of any gate effect, the effective density of localized states available for hopping transport 
at some temperature T is given by z/(e ~ kBT) where energy is measured from the band 
edge. An exponential DOLS of the Vissenberg model would then lead to an exponential 
dependence of /io on T. Gate voltage dependence in this case comes from electrostatic 
modulation of the Fermi level in this rapidly varying DOLS. It is difficult to understand 
otherwise how an exponential dependence of the prefactor on T or Vg could arise. Deriving 
a quantitative relationship between Vg and the local Fermi level would require solving the 
full electrostatic problem of SCPLF conduction in the presence of the transverse gate field. 

One can consider whether the proposed SCLPF conduction is a bulk process or one 
dominated by conduction in the thin channel layer active in standard FET operation. For 
the SCLPF mobility parameter fiQ to coincide with the zero field mobility found at high 
temperatures in the linear regime, an effective cross-sectional area A considerably larger 
than wx a few nanometers is required. This is also true for the data of Ref. j2^, in which 
a considerably different device geometry was used, if analyzed using the SCLPF model. 
However, the model of Eqs. (H)),© does not account for the presence of a gate electrode or 
complicated source and drain geometries, and full computational modeling in this regime 
may be required to quantitatively account for this. One test for bulk vs. channel conduction 
would be to search for a correlation between P3IIT thickness and currents in the apparent 
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SCL regime. This investigation is ongoing. 



V. CONCLUSIONS 

For a series of field effect devices witfi cfiannel lengths ranging from 3 /xm to 200 nm, we 
find gate modulated nonlinear conduction at low temperatures and high average source-drain 
electric fields. Analysis of channel and contact resistances as a function of temperature, and 
the dependence of conduction on sample size at low temperatures support the conclusion 
that this nonhnearity is unlikely to be a contact effect. We find that the data are well 
described by a model of gate modulated space-charge hmited currents with Poole- Prenkel-hke 
behavior of mobihty. Within this model, the Vq and T dependence of the mobihty prefactor 
is consistent with a very strongly energy dependent density of localized states. Finally, 
the temperature dependence of the Poole-Frenkel-like term within this model suggests a 
crossover from thermal hopping to quantum tunneling at low temperatures. Further studies 
of the field effect electrostatics problem, the metal-semiconductor contacts, and the low 
temperature nonthermal hopping process should lead to increased understanding of the 
conduction processes at work in these materials. 
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FIG. 1: Mobility vs. T as computed from transconductance for three gate voltages, for a device 
with L = 300 nm from room temperature down to 200 K. Lines are fits to various models of 
hopping transport described in the text. Inset: cross-section of device showing bottom contact 
configuration and definition of channel length, L. 



FIG. 2: Id vs. Vd for the L = 3 fim, w = 30 fim device, at 300 K (nearly linear, top) and 70 K 
(highly nonlinear, bottom). Curves from the top down correspond to Vq values from -95 V to 
-30 V in intervals of 5 V. 



FIG. 3: (a) Log-Log plot of Id vs. Vd for a device with L = 500 nm at Vq =-75 V. Solid lines 
indicate a numerical solution assuming space-charge limited conduction with a Poole-Frenkel-like 
field dependence of the mobility, (b) Plot of ln(//5/V^ vs. \/Vd, as suggested by Eq. 13 Solid lines 
are fits to a linear dependence on \/Vd- 



FIG. 4: Plot of Ron X w vs. L for the fixed-aspect-ratio device series in the linear regime at 300 K 
for several gate voltages. Slopes of the linear fits correspond to i?^ of the channel, while intercepts 
correspond to the parasitic contact resistivity, RgW. Inset: the ratio Rg/Rch for the L = 1 //m, 
It; = 10 iim device. Clearly this device is not contact limited at room temperature. 



FIG. 5: Rs/Rch as a function of temperature as determined for the w = 100 /xm, L = 5 /im member 
of the fixed- width set of devices. Since this ratio decreases as T is lowered, contacts actually improve 



relative to t 
elsewhere 



le channel at low temperatures. This interesting result will be discussed more fully 



24|. 
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FIG. 6: Inset: Plot of the parameter 7 vs. Vq for the 500 nm sample of Fig. El at various 
temperatures ((top-to-bottom) 10 K, 30 K, 50 K, 70 K, 90 K, 120 K, 150 K, 180 K, 210 K) 
showing that 7 is roughly gate voltage independent. Main figure: Plot of 7 vs. 1/T for several 
samples, with 7 averaged over gate voltages for each sample. Error bars are standard deviation. 
At high temperatures 7 is expected to vary linearly in 1/T, and all samples show a similar slope, 
7 X T ~ 0.12 (m/V)^/^. Within this model, saturation of 7 at low temperatures would indicate a 
crossover from thermal to nonthermal hopping transport at low temperatures. 



FIG. 7: Assuming a fixed effective area A = 3 x 10^^^ m^, (a) Plot of the parameter hq from plots 
like Fig. vs. Vq for the 300 nm sample, from 210 K down to 10 K. (b) Plot of fiQ vs. T for 
the same sample, for 5 different gate voltages. Note that //q depends nearly exponentially on both 
temperature and gate voltage. 
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